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abstract The increasing rearrangement of a univariate function is the increas- 
ing function which yields the same image of the Lebesgue measure than the orig- 
inal one. The convex rearrangement of a function is obtained by integrating the 
monotone rearrangement of its derivative. Both operators can be generalised 
to higher dimensions, where an increasing function is the gradient of a convex 
function. We define here the rearrangement of an irregular function as the limit 
of rearrangements of approximations, and give a consistency theorem. 

In this paper, we investigate the asymptotic rearrangements for approxima- 
tion of random fields. Stronger results are given for Gaussian fields, and the 
examples of Levy fields and Chcntsov fields are derived. 

keyword random fields, rearrangement, limit theorems, random measures 

MSC 60G60, 60B12, 60G57 

Introduction and notation 

The following notation will be useful. In R d , denote by e = (e l )i<j<<j the 
canonical basis, and + stands for the Minkowski addition of sets. The operators 
vol, diam, cl, int, d resp. stand for the volume, diameter, closure, interior and 
border of a Borel set. For a centred random vector Y taking values in K d , cov(Y ) 
is its covariance matrix, with cov(Y)ij = M(YiYj), 1 < i, j < d. Given two real 
random variables Xi, X2, also write cov(Xi, X2) = E(XiX2) — E(Xi)E(X2) their 
covariance. The weak convergence of measures is denoted by =>. 



Preliminary example Consider a finite population, arbitrary labelled with 
numbers k in {1, N}, for N G K* . For 1 < k < N, the member k receives 
an income of a certain resource, denoted by a real number g(k). Now call a 
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the permutation of {1, N} which makes the function k — >• g(a(k)) increasing. 
Call j=joit the monotone rearrangement of g. 

Define C(k) = 2i=iff(*)> 1 < k < n. Since g is monotone, C is convex. 
For 1 < k < N, G(k) represents the total amount of resources detained by the 
— -th poorest fraction of the population. Now, call C(k) = —C(n). It is the 
"equality function" , in the sense that C = C iff all incomes are equals. Also, for 
some distance 5, the distance Is = S(C, C) between C and its equality function 
measures the inequalities among the population. 

If one defines f(k) = 2»=i5(*)j 1 < k < N, the cumulated income, C is 
called the convex rearrangement of /. It is indeed the only convex functions 
which has the same increments (but in an other order), and coincides with / in 
1 and N. Consider for instance the case where 5 is the L 1 norm on M. N . For a 
given cumulate income function /, the quantity 

N , 

Is = £(C(fc) - -f(n)) 
k=l 

retrieves the Gini coefficient, which has played a central role in measuring eco- 
nomic inequality since its introduction by Corrado Gini at the beginning of the 
20th century. The use of the convex rearrangement for measuring economic 
inequality is discussed in [12]. 

The notion of rearrangement, defined above for a discrete population, can 
be generalised in the continuous framework. If g\ is an integrablc function on 
[0, 1], and a is a transformation of [0, 1] which preserves Lebesgue measure, the 
function defined by 

92 = 9i o a (1) 

is a rearrangement of g. Roughly, the function g^ takes each value the same 
number of times than g\ does, but not at the same places. For any function g, 
denote by [i g the image of Lebesgue measure under g. Relation (fT]) also implies 

M S2 = /%i ■ ( 2 ) 

A function g2 is said to be a rearrangement of g\ if it satisfies ([2]), but this is 
not equivalent to (JlJ. Consider for instance g\{x) = x and g2{x) = (2x — [2a:]) 
on [0, 1]. They are rearrangements of one another but it is not possible to find 
a such that ([T]) is satisfied. 

The situation on a compact body K of M. d is similar. For a function g in- 
tegrablc on K, call /i g the image of Lebesgue measure Xd under g. Then a 
function 52 is a rearrangement of an other function g\ iff it satisfies ([2]). The 
rearrangement is furthermore said to be monotone if 172 is a monotone function, 
i.e the gradient of a convex function. According to Brenier's theorem, stated in 
section [T] each integrablc function admits a unique monotone rearrangement. 
The result includes a non-degeneracy conditions on \i g for having ([T]) as well. 
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In this article, we define the convex rearrangement of a multivariate smooth 
function / as the function obtained in rearranging the gradient of / monotoni- 
cally. It generalises the notion of convex rearrangement on {1, N} introduced 
in the preliminary example. 

If a function is irregular, one can take regularisations and study asymp- 
totically their rearrangements, under the proper renormalisation. One of our 
main results states that the asymptotic rearrangement is consistent, i.c if fi 9n 
admits a weak limit, for a sequence of functions {g n ; neN}, then the mono- 
tone rearrangements of the g n also converge. The result is similar concerning 
convex rearrangements. With this tool, asymptotic convex rearrangement of 
some random fields are studied, with stronger results for Gaussian fields. A 
section is devoted to the study of examples, like the Levy and Chentsov field. 
It will become apparent in this paper that asymptotic rearrangement is a com- 
plex transformation for multivariate fields. In particular, the results are highly 
dependant on the method of approximation. 

It is of practical and theoretical interest to investigate asymptotic proper- 
ties of rearrangements. It can be used, for example, to construct estimators 
of parameters of stochastic processes, and for measuring their fluctuations, see 
[7J. There are also connections between convex rearrangement and other areas 
of research such as Finance Mathematics and Economics. The Lorenz curve, 
important in finance mathematics, is a common object in convex rearrangement 
of Gaussian processes. In the field of econometrics, convex rearrangement can 
be used to measure the indices of fluctuations of stochastic processes, related 
to indices of economic inequality , like the Gini index in the preliminary exam- 
ple, see [12j . The monotone rearrangement of a function g also has a physical 
meaning, as the solution of the optimal transport problem with transfer plan g. 

The first section is devoted to explain the principles of rearrangement in one 
dimension. The example of Brownian motion is developed as explanatory exam- 
ple. A survey of the results concerning rearrangements of stochastic processes 
of one variable can be found in [7J. Then, a reduced version of the problem 
of optimal transport is introduced, and we give its connection with monotone 
rearrangement. An adapted version of Brenier's theorem is given, which yields a 
rigourous definition of monotone and convex rearrangements. A necessary and 
sufficient condition for the convergence of the rearrangements of a sequence of 
functions is the main result of this section. 

In Section[2j given a family of random fields {Y n ; n <G N}, typically irregular, 
we investigate conditions that ensure the existence of the limit of their rearrange- 
ments. We obtain stronger results in the case of Gaussian fields approximated 
by polygonal interpolations. The efficiency of the method is illustrated in Sec- 
tionO where various types of examples are derived using the convergence results. 
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1 Monotone rearrangements and optimal trans- 
port 

This section exposes the theoretical material required for rearranging multi- 
variate functions with compact support. It is related to the optimal transport 
problem, in that the monotone rearrangement 9JtT of a transport plan T coin- 
cides with the optimal solution to the corresponding transport problem. Then, 
we study the consistency of the monotone rearrangement, needed for rearrang- 
ing irregular functions, the same way it is done for Brownian motion just below. 



1.1 One dimensional case. Convex rearrangement of the 
Brownian motion 

The monotone rearrangement of an intcgrablc function 5 on [0, 1] has been de- 
fined in the introduction. It is the unique increasing function with the same 
distribution, and is denoted by 9Hg. It has been shown in [§] that every in- 
tcgrablc function on [0, 1] admits a monotone rearrangement, unique up to a 
negligible set. We emphasise here that the central object of the monotone rear- 
rangement is /Lt g , the image of Lebesgue measure under g. In other words, two 
functions have the same rearrangement if they have the same distributions. 

If now / is an absolutely continuous function on [0, 1], i.e such that for all 
x in [0, 1], f(x) — /(0) = J g(t)dt for some integrable function g, the convex 
rearrangement of / is the unique convex function C verifying C(0) = /(0) and 
C = a.c. Write C = £/, where £ is the convex rearrangement operator. 

For / irregular, one chooses smooth approximations {/„ ; n > 1}, and 
studies the asymptotic of the rearrangements. If there exists a sequence {b n ; n > 
1} and a convex function C such that -^-<tf n — > C a.e, C is said to be an 
asymptotic convex rearrangement of / with renormalising sequence {b n ;n > 1}. 

There is a wide literature investigating the asymptotic convex rearrange- 
ments of random processes. Although a rigourous study is not trivial, it is 
possible to understand better the convex rearrangement machinery in the case 
of the Wiener process. Take X a standard Brownian motion on [0,1], with 
X n its piece- wise linear interpolation on < fc < n}, normalised by y/ri to 
avoid the divergence of its derivative. For each n, X n is differentiable a.e., and 
the image of Lebesgue measure Ai under the renormalised derivative is writ- 



ten /!„ = Ai I -j=X' n ) . The independence of increments implies that /j, n is 



the empirical measure of n independent normal variables, and it is clear that 
it will converge weakly to the normal distribution 71. It will be rigourously 
proven later, in theorem II .21 why this implies that the asymptotic convex rear- 
rangement of X on ]0, 1[ is the Lorenz curve GL\, defined as the unique convex 
function with gradient distribution 71. Davydov and Vershik [B] obtained the 
strongest result, namely the uniform convergence of \\^€X n — GLiW^ to 
with probability 1. 
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A lot of similar results are obtained with processes that have stationary in- 
crements, or are stable, see the survey [7]. Azais and Wschebor Q] also showed 
that, for AT in a certain class of Gaussian processes, if instead of a piece-wise 
linear approximation, one chooses for X n a regularisation of A by a convolu- 
tion kernel, then A admits the same asymptotic convex rearrangement, namely 
the generalised Lorenz curve GL\. In this case, the asymptotic convex rear- 
rangement of / seems unambiguous, modulo the multiplication by a non-zero 
constant, in the sense that it does not depend on the approximation method. 
We will see in section[3]that it is not the case for most multivariate random fields. 



1.2 The optimal transport problem and rearrangement 
operators 

The problem described below is a simplified version of the traditional optimal 
transport problem, which is fully described and exhaustively discussed in jTTj . 

A company has a capacity of production per unit time represented by a 
measure fj, on M. d , the production measure. The quantity produced in area dx 
per unit time is [i{dx). This company has to deliver its products to a domain K 
of K d , compact and convex, where the demand is uniformly distributed. The cost 
of transport between a site of production s and a point z in K is proportional 
to || s — 2 1|, the Euclidean distance. A transport plan S is a function which 
associates to each z in K the corresponding production site S(z), where the 
product delivered to z comes from. In particular, we need to have, for all Borel 
set B <= B d , 

H S (B) = KB), (3) 

so that all production sites provide the required quantity of product. The total 
cost of this transport plan (per unit time) is hence 

C(S) = [ \\z - S(z)\\dz. 

JztK 

The optimal transport problem consists in finding a transport plan S : K — > R d 
minimising the cost C(S) under requirement Addressing this issue, suppose 
that a given transport plan S is modified by switching the destinations z and £ 
for two productions sites S(z) and £(£) for an infinitesimal quantity of product. 
The new transport plan is denoted S and the corresponding cost variation is 

C(S) - C{S) = 2(z - C, S(z) - S(())(dz + d(). 

Informally, a transport plan will be in some sense locally optimal if, for all 

(z-C,S(z)-S(O)>0. (4) 
It turns out that (J4j) and <j3j> indeed characterise optimal transport plans (see 

HQ)- 
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The question that naturally arises now is about the existence of such an 
optimal transport plan. That is the purpose of the following theorem. 

Theorem 1.1 (Brenier). 

Call IC(K) the class of convex functions on K . LetQ(K) be the set of monotone 
functions on K , defined by 



Q{K) = jvC | C G K(K), J || VC|| < ooj . 

Then, if [i is a measure on R d with finite first moment, there is a unique mono- 
tone function in Q{K), denoted by M^, such that A^M" 1 = a. If a point zq 
of K is unambiguously defined as "starting point", call C M the convex function 
which gradient is M^, satisfying C^(zq) = 0. 

Comments, proof, and a more general result can be found in [3]. The function 
M M is the optimal solution of the transport problem with production measure 

Theorem 11.11 is the proper tool to define high dimensional monotone and 
convex rearrangements. 

Definition 1.1. For an integrable function g on K, define SSJlg = M t , s its 
monotone rearrangement. 

For a function f which gradient is integrable, there exists a unique convex func- 
tion, denoted €f, which satisfies 



X d (V€f)- 1 = A d (V/)-\ 

e/fo) = /(*>)• 

It is called convex rearrangement of f. 

The convex rearrangement can also be defined as £/ = C pv/ + f(z ). Given 
a vector- valued function S on K , since VJIS is the gradient of a convex function, 
its restriction to each segment [z,(] c K is increasing, and hence it satisfies (|4]). 
In this regard, theorem 11.11 provides with a unique solution to the optimal 
transport problem with transport plan S. Note that [3] also gives the exis- 
tence of a measure-preserving transformation a of [0, 1] such that QHS* o a = S ' , 
provided (1$ is absolutely continuous with respect to Lcbcsgue measure, which 
justifies the "rearrangement" terminology 

In dimension 1, convex rearrangement was already defined in the literature. 
The class S(K) is exactly that of absolutely continuous functions if if is a 
compact interval of R. Hence, it is a generalisation of absolutely continuous 
functions upon which we extend operator £. Note that, although it is called 
"convex rearrangement" , function £/ is not a rearrangement of / in the sense 
of P]). For instance, / and <£/ do not in general yield the same maximum. 
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Nevertheless, visually it corresponds in some way to piling up the increments of 
/ in another order. 

Here is a simple example that allows us to rearrange multivariate hclds with 
a simple structure. 

Proposition 1.1. Let A\, ...,A d be compact intervals o/M and Jor each i in 
{1, d}, let ft be a function in S(Ai). Then, the function f = © i=1 fi, defined 
on K = Ai x ... x A d by 

d 

Qf i (z) = f 1 (z 1 ) + Mz 2 ) + ... + f d (z d ) 
1=1 

satisfies 

feS{K), 
a 

e/ = 0<£/i. 

i=i 

Concerning the starting point zq, it is arbitrary and plays no role. It will 
sometimes be implicitly defined and not mentioned, but one has to be consistent 
with its use. For instance, in the previous example, operator £ is defined on 
subsets Ai of K, with implicit starting points zo.i. In order for the assertion 
to be true, one has to take on K as starting point of convex rearrangement 

z = {Z0,i)l<i<d- 

Proof or vrovosition \l.l\ Denote by (^) the product of measures. One simply 
has to use the fact that 

d d d 

mv/ = = ® = 0M(c/i)'> 

i—l i—1 i—1 

since the SOT /■ are rearrangements of the fl,i = l,...,d. We also have, 

d 

i=i 

and this means that ©i=i i s a convex function having the same gradient 
distribution as V/. It yields the result because these functions coincide in zq. 

□ 

1.3 Consistency of the rearrangement operators 

In this article we deal with irregular random fields, for which we cannot a priori 
obtain a convex rearrangement due to the absence of gradient. In consequence, 
by analogy with the 1-dimensional case, we instead investigate asymptotically 
the convex rearrangement of their regularisations. Call asymptotic monotone 
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rearrangement of g any monotone function that is the limit of renormalised 
monotone rearrangement of smooth approximations of g. The theorem II .21 will 
allow us to obtain an asymptotic monotone rearrangement of a function study- 
ing the asymptotic of its gradient distributions. 

Similarly for a function /, call asymptotic convex rearrangement of / the 
limit of convex rearrangements, after renormalisation, of approximations of /. 
It will become apparent along this article that, if one simply uses this definition, 
the set of all possible asymptotic monotone rearrangement of a given function 
is wide. To have consistency between convex rearrangement and asymptotic 
convex rearrangement for a smooth function, one need at least to control the 
gradient of the approximation. 

Denote by u n A-it the point-wise convergence on the interior of K. If K is 
implicit, it will be omitted. In the sequel, if is a convex body of M. d , with an 
arbitrary starting point zq G K . The following theorem will be our main tool 
for rearranging random fields. 

Theorem 1.2. 

Take {/„; n > 1} and f in S(K), and define g n = V/„, g = V/. Then the 
three following statements are equivalent: 

M ff „ %, (5) 

Mg n -> Wig, (6) 

e/n-/»(*b) -> £/-/(%). (7) 

The proof is in section IQ1 

Remark 1.1. The convergence here only occurs on int(K). In order to obtain 
uniform convergence on all K , one has to control that the gradient is never "too 
large" on "small sets", typically located on the edges of K . Formally, we have the 
following sufficient condition, proved in JSjj. With the previous notation, note the 
gradient's norm distribution function by F n (a) = X c i({z; \\V f n (z)\\ > a}), a > 0. 
If 

(Fn) 1 ^ is uniformly integrable in + oo, (8) 

then (£/n) is equicontinuous on K , and the convergence is uniform on K . Due 
to the light tail of Gaussian distributions, it is likely that the Gaussian fields 
under study will satisfy (Sj) after approximation and renormalisation, but it is 
not the purpose of this article to obtain optimal results, and we only checked 
condition |5|] in our examples. 

A similar result was already proved in the 1-dimensional case in [J , including 
the corresponding condition for uniform convergence on all K, not only its 
interior. In the probabilistic framework, in order to check condition ([5]) for 
random fields, it is easier to deal with characteristic functions. We use the 
following lemma: 
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Lemma 1.1. Let (fi n )n>i be a sequence of random probability measures with 
characteristic functions (tp n )n>i- If there is a probability measure /i with char- 
acteristic function <p such that, for all h in M. d , 

fn(h) — > ip(h) a.s., 

then \i n =>■ fi with probability one. 

Proof. We have 



^{ t p n (h,u)-Kp(h,u)}dh¥(dw) = 1. 
Due to Fubini's theorem, with probability one, for almost all h of M. d , 

Ifin(h) -> ipQl), 

and it is well known that it implies the weak convergence of the corresponding 
probability measures. □ 

The next result can also be used. Following [2], call convergence-determining 
class C a class of Borel sets such that the weak convergence of measures follows 
from the point-wise convergence on C. Theorem 2.2 p. 15 in [Sj implies that 
there is a countable such class in R d , and hence in our case it is enough to show 
Hn{E) —> fJ,(B) a.s, for all /i-continuous Borel set B in the Borel cr-algebra Bd 
of R d . 



2 Asymptotic rearrangement of random fields 

In this section, we consider a random field X defined on Kd and give general 
results about its asymptotic rearrangement. Then we give the main theorem 
of convergence in the case of Gaussian fields, in the framework of polygonal 
approximation. This generalises the asymptotic convex rearrangement of the 
Brownian motion derived in Section ll.il 

2.1 General results 

The notation {Y n ; n > 1} stands here for a sequence of smooth vector valued 
random functions, and {/i„ = fj,y n ] n > 1} are their distributions. In this 
section general results concerning the asymptotic of (/i n ) n >i are given. The 
objective is to obtain a deterministic limit measure [i of the /i„ , and use the 
consistency theorem 11.21 We hrst derive an expression of the only possible limit 
li with the help of Fubini's theorem: Given a measurable set B of Bd, we have 

E(> n (B)) = E ( / l Yn{z)eB dz) = f ¥(Y n (z) g B)dz. (9) 

\JK d / JK d 

Letting n go to oo in © yields the following: 
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Proposition 2.1. We assume that the two following conditions hold: 

There exists fi probability measure on W 1 such that, a.s, pL n =>■ fj,, (10) 

For a.e. z in Kd, 

there exists [i z probability measure on R d such that fJ.y n (z) =^ l^z- (11) 
Then, 

fx(B) = [ n z {B)dz, BeB d . (12) 

We also have ?fflY n — > M M in virtue of theorem \1.°A 

Now, as a first example, the following proposition gives a sufficient condition 
on the conjoint laws of the variables (Y n (z)) z eK d for the convergence of jj, n . 

Theorem 2.1. Suppose that ill]) is satisfied. Define /i by llfy) . Let Y a i n \ be a 
sub-sequence such that, for all /z- continuity sets B in a convergence-determining 
class of Bd (see. JBjj, p. 15), 




then /x CT („) M a - s - 

Wc hence have asymptotic convex rearrangement with probability 1 for some 
sub-sequences. 

Proof. Without loss of generality, we assume that cr(n) = n. For B a /i- 
continuity Borel set in the convergence determining class, 

E (\Hn(B) - E( Pn (B))\ 2 ) = E {^ n (B) 2 ) - (E Mn (S)) 2 
= e( [ dzl Yn{z)eB [ d(l Yn(0eB )-[ dz¥(Y n (z)eB) [ d(P(Y n (0eB) 

\JK d JK d J JK d JK d 




Hence, hypothesis ([TU)) , along with Borel-Cantclli's lemma ensures that for 
all B in Bd, with probability one, ^i a ^(B) — > fJ-(B). This relation is hence 
true a.s simultaneously for all Borel sets of a countable generating subclass of 
Bd, and so fJ, a ( n ) /i with probability l(See theorem 2.2 in [2]). Theorem 11.21 
completes the proof. 

□ 

For most of the random fields investigated in Section[31 cov (ly„( z )g_B, ^-Y n (QeB) 
is in O(i) and we cannot have asymptotic rearrangement for 9}lY n . We need 
stronger results in this case, and were able to obtain them in the framework of 
Gaussian fields, interpolated on a simplicial triangulation. 
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2.2 Simplicial approximations on Kd 

Most of the commonly investigated random fields of the literature are irregular, 
and hence cannot be directly rearranged, they need to be approximated by 
smooth functions. In this article, we only followed the following paradigm: 
Given a random real field X, define approximations X n of X, then normalise 
and rearrange monotonically their gradient, which will be called Y n = r-VX re . 

In this paradigm, one would like the result not to depend on the choice of 
the approximation X n , as long as it converges to X. It is in the very nature of 
the convex rearrangement to be sensitive to slight changes in the approximation 
method. Consider for instance the following deterministic example. Define /„ 
as the continuous function on [0, 1] null in 0, linear on each segment [— , ls±L^ f or 
1 < k < n, and with slope ±1. Then, /„ uniformly converges to the (convex) null 
function, but <£/„ uniformly converges to the convex piece-wise linear function 
null in having slope —1 on [0,^] and +1 on [|,1]. To avoid this kind of 
phenomenon for asymptotic convex rearrangement, one needs to control that the 
gradient of the approximation resembles the gradient of the original function, 
or its variations if there is no gradient. That is one of the reasons why we chose 
for X n the interpolations of X on a reasonable subset, namely the vertices of a 
triangulation. We present now the definitions and notations needed for stating 
the results. 

Denote by Kd the unit cube [0,l] d . Call simplex of W 1 the convex hull 
of any (d + l)-tuple of points with non-empty interior, and call triangulation 
of Kd any finite simplicial partition of Kd- Let T be such a triangulation 
of Kd- Denote by X T the simplicial approximation of X with respect to T, 
i.e. the function which is affine above each T in T and coincides with X 
above the vertices of T . If (7^) n >i is a sequence of triangulations such that 
diam(7^) = sup T67 - n diam(T) — > 0, it is called an approximating triangulation. 
In this case denote by Xj the corresponding approximation, i.e X% = X Tn . We 
will consider in this paper exclusively approximating triangulations of a special 
form, described below. 

Call germ of triangulation any finite set of simplexes T verifying the following 
property: There exists a network T of M d such that 



Any network T satisfying (|14p is said admissible for T, and the notation IV 
refers to an arbitrary choice of such a network. Then define, for n > 1, 



Property (fl"4")) ensures that T n is indeed a partition of Kd. The problem is that 
the set — (•j + T) n Kd might not be a simplex. However, those problematic 
— (7 + T) won't play any role in the asymptotic convex rearrangement because 
their number is negligible. So, arbitrarily divide each of them in a simplicial 



{7 + T I 7 G T, T e T} is a partition of R d . 



(14) 
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partition such that there is, for all n^a, triangulation T n which is a simplicial 
sub-partition of T n , and differs from T n only regarding the simplexes touching 
the border of Kd- Identify T with the simplicial approximating triangulation 
sequence (T n ) n >i- 

Since X% is a.e afhne, denote by VX^~ its gradient. In all the paper, 
{b n 'i n > l} n >i stands for a sequence of positive numbers which aims to give 
sense to lim„ j-9UVXj" (or, cquivalently- see theorem II .2[1 - to lim n ^-<ZX[[). 
The renormalised gradient is denoted by 



Y T = 



VX 



T 



Using theorem II .21 to obtain the rearrangement of Y^f, it is more convenient to 
work with its distribution = /jyr . 

Working with random fields, given a germ of triangulation T, it is more 
feasible to study the behaviour of X^ above the translates of only one simplex 
at a time. We state in this purpose the following lemma: 

Lemma 2.1. Let T be a simplex ofT . Put 



r (T) 

HP 

,,( t ) 



-h + T)c K d 
n 



T + lt T > 



Y n (z)l {zeH (T)y z G K d , 



Then, for all B in Bd, 



TET 



Proof. For r > 0, call B(0, r) the open ball with radius r and centre 0. 



TeT 



< 



< 



12 1 {±h+T)ndK d m Xd (^ +T n 



TET,iET T 



< X d [dKd + --8(0, sup diam(T) ] 



A 
C r 



TET 



(15) 



where Cf is a constant depending on T. To be convinced of the last formula, 
one needs to notice that it is possible to find q n = Cdn d ~ x points on dKd 
z 1 , z qn (for some constant Cd) such that OK + ^-E>(0, Cr) is always included 
in U£i B(z l , 2C-r/n), where + is the Minkowski addition of sets. □ 
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In conclusion, for the study of a random field X with a germ of triangulation 
T, the proceeding is the following. For each T G T, compute the measure /A 
by studying separately the increment along each edge of simplex T . If the limit 
// T ) exists for each T, then the limit rearrangement is given by the theorem 1 1.2 1 
with limit measure = X^tgt • 



2.3 Rearrangements of centred Gaussian fields 

The specific study of Gaussian fields yields more efficient tools to study the 
convergence. We give here the statement of the main theorem of this section, 
some examples will be derived in the next section to illustrate the theory, for 
multivariate Levy and Chentsov fields. The generalised Lorenz curve plays a 
great role in the convex rearrangement of Gaussian processes, so we introduce 
it now. 

Definition 2.1. Call 7^ the d- dimensional standard normal distribution. The 
d-dimensional generalised Lorenz curve is 



GL d — C ld . 

In other words, it is the asymptotic convex rearrangement of any field which 
renormalised gradient measure converges to 7^. It corresponds in dimension 
1 to the classical Lorenz curve, frequently used in the fields of finance and 
econometrics. 

Following lemma |2~T1 choose T a simplex of T. Here Y 7l is the gradient of 
Xn . Call A„ (z) — cov(y„ (z)) Let us translate the results of proposition 
I2.1l for Gaussian fields. 

Proposition 2.2. Assume that \10\) is fulfilled and that for all z there exists a 
matrix hS T > (z) such that 

A^(z)~,A^(z) (16) 

(which is equivalent to Ul\) in the Gaussian case). Let y} z be the Gaussian 
measure onW 1, with covariance matrix hS T '{z),z G K c i- Then, the only possible 
limit probability measure pS T ^ is the Gaussian mixture 

^ T \B) = [ lxf\B)dz, BeB d . 

We state now the main theorem of this paper, which gives a more efficient 
condition for the convergence of p^h"^ = /i y <T) than theorem 12. II We first state 
the technical optimal result, and then give the corollary that will be our main 
tool for rearranging random fields. 
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Theorem 2.2. Let X be a Gaussian centred field on Kd, n € N, and T a 

regular simplex of Kd admitting u = (u i )^_ 1 as orthonormal basis. Define Hn 

(T) (T) 

as the distribution of Yn and ip n is its characteristic function. Then there 
is a deterministic sequence c n of order O (i) that satisfies 



E(\^(h)-E(^(h))-c n \ 4 ) 
= o([ dzdCe^^ A " T, ^ +A " T> ^ sup \E(Y n Az)Y n ,(C))\] ,heR d . 

\J{K d )2 l<i.j<d J 

(17) 

The proof is in Section HT21 In all the cases treated here, this result gives us 
the convergence of the sequence of random measure fj,^ to the limit if it 

exists. The deterministic term c n is an effect of the problems occurring at the 
border, see lemma I2TT1 

A more handful version of previous result is the following corollary, which 
gives the convergence for many random centred Gaussian fields with piece-wise 
smooth covariance function. 

Corollary 2.1. Let X be a centred Gaussian field, and T a germ of triangu- 
lation. Call T the covariance function of field X. We set the set of all pair 
of points in (Kd) 2 upon which T is not of class C 1 . If U6}) is fulfilled for each 
simplex T in T and the following holds, 

n « 

E(^l eni ^|) (19) 
sup ||Al T )(z)|| <oo, (20) 

then we have, with probability 1, 

fil^fi T , 0JtY n —> M M r, ^<LX n -> C M r. (21) 

On 

where fj, T = X^tsT A* ■ 

Proof. Take T e T, and u a basis of T. For (z, () 9, call dr z ^ the differential 
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form of r in (z, Q. We have, 

Vl<i,j<d, (z,C)e(ffi T) ) 2 \e, 



6 2 



1 



1 



n n 
1 , 1 



1 



T[z + -u\ c + -u> - r * + -u l , c -r Z) ( + -rf +r(z,() 



l 



SK c U ui '^)- dr *- c u ui ' 0) J- dr *4°^ u3 



= o 



The others (z, C), where T is not regular, arc supposed to be negligible, and we 
have the following upper bound : 

which is uniformly bounded by hypothesis. In consequence, due to the previous 
theorem, with t n = vol (^T 1 ), 



(h,(cov(Y^ T » (z))+cov(Y^"> (C)))h) 



sup 
i-j 



( \ 



dzdC 



We easily check that t\ — O (-^j) and so, using theorem 12. 2\ we have, for a 
constant C 



1 



,2il 



\ 



< CO. 



7 



Due to Borel-Cantelli's lemma, a.s (/?„ '(h) — E(</?„ (/i)) — c„ — ?> 0, for each 
T e T- (Remember that c n is a deterministic sequence in 0(l/n)). Since by 
hypothesis E(^ T) (h)) -> (^ (T) (/i) for each T, we have the convergence of </3„ (/i) 
to (p(h) with probability 1, for each h in R d . Lemma [Ol and theorem 1 1 . 2 1 carry 
the conclusion. 

□ 
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Theorem 12.11 only allowed an upper bound of order i, but with the last 
expression one is often able to derive a summable upper bound. The scope of 
this article is the study of Gaussian fields, but one could probably compute the 
quantity E (\<p n (h) — E(</?„(/i))| 4 ) for a wider class of fields and obtain an upper 
bound better than — , at least for processes with low-dependant increments. 



3 Examples 

Write S d = {(U)i<i< d \0<U< l,£ti U < l} the elementary simplex of M . 

Given z in M. d , an orthonormal basis u = (u l )x<i<d of M. d and a positive number 
I, define the regular simplex with summit z, basis u, and side-length I as 

T(z,u,l) = z + lp u (X d ), 

where p u is a rotation of M. d transforming e into u. If u = e, T is said to be 
"straight". If I = 1, T is said to be "regular". 

There is a large literature in the 1-dimensional case. Corollary 12.11 enables 
us to retrieve partially some results. 

Theorem 3.1. Let X be a 1-dimensional Gaussian centred process with sta- 
tionary increments, put a 2 (t) = E(X(t) 2 ). We set b n = na(-) and make the 
following assumptions 



a 2 



4 < oo, (22) 

n n 

is piece-wise of class C . (23) 



The 



^€X n ^GL 1 . 

On 

Proof. By stationarity 

E(i:(z) 2 )=E(I,',(0) 2 ) = «V 



and so we put Y n = ^-X' n . Then, for all z, cov(Y n (z)) = cov(F„(0)) converges 
to 1, and (fTB")) is satisfied with A(z) = 1. So, the candidate for the limit, given 
by (fT2|) . is 71. Also, sup„ gN |covy„(0)| < oo and ([20]) is satisfied. Condition 
(fT5|) is also satisfied by hypothesis. To check (|19p . we call E the finite set of 
points upon which a 2 is not of class C 1 . We have 

e = {(s,s + y) | s€[0,l],y€E}. 

Now, n H 2 \ = \E\ \H n \ = 0(n), so (fT9| is also satisfied, and a.s yU„ => 71. 
Now, GLi is indeed the convex function whose gradient distribution is 71, and 
theorem 11.21 gives the conclusion. □ 
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Davydov and Thilly [5], along with a general theorem concerning Gaussian 
processes with stationary increments, have obtained the following result for the 
fractional Brownian motion: 

Theorem 3.2 (Davydov, Thilly 98). Let < a < 2, and W a be the standard 
fractional Brownian motion. Then 

\\-^<LW:-GL4W^0. 

In this case, the covariance function cr 2 (t) = \t\ a is non differentiable only 
in 0. We set b„ = ?i 1 ~ Q / 2 , and is summablc iff a e]0, §[. Hence theorem |3. II 
retrieves the result for < a < 3/2. This gives us the convex rearrangement 
of Brownian motion, first stated by Davydov and Vershik [5] , who furthermore 
obtained uniform convergence on [0,1]. 

Now we present some multi-dimensional examples, which give an idea of the 
variety of possible phenomena that can occur. 



3.1 Levy field 



The Levy field is the first multivariate random field that we study. In dimension 
d > 1, it is defined as the only centred Gaussian random field with covariance 
function 

r(*,c) = 5(NI + ||CI|-|k-CII). 

It is spherically symmetric and its 1-dimensional version is the standard Brow- 
nian motion. 

Let X be a Levy field. We use the notation of section 12.21 Let us take 
an arbitrary simplex T = T(zq, u, 1) in (R + ) d . Then, the covariance matrix of 



VXn \z) is given by 




+ 0(1)) 



^ + o(l)) ifi^j, 



A first observation is that the covariance matrix does not depend on the choice 
of orthonormal basis u. It is not surprising, given the spherical symmetry of 
the field. But since the covariance matrix was computed in basis u, the limit 
distribution does depend on u. We set 



C=1 _V2 
2 ' 



A = 



/ 1 

C 

V o 



c 



c 



c \ 

c 

1 / 



17 



and let p be the Gaussian distribution with covariance matrix A in basis u. 
The matrix A admits 1 — C as eigenvalue of multiplicity d — 1, associated to 
cigcnspacc (1, 1) , and [n — 1)C+1 as eigenvalue of multiplicity 1, associated 
to eigenvector (!,...,!). 



Since 



e n {Hk 



0(n d ), using corollary 12. II and theorem 11.21 we have 



the three usual results 



The asymptotic rearrangement is consistent under the action of rotations: In- 
deed, if pT is the limit measure with germ of triangulation T, we have for all 
rotation p and germ of triangulation T, This is due to the 

symmetry property of the Levy field, and will not be the case in the subsequent 
examples. 



3.2 An additive Gaussian field 

Define X(x, y) = W(x) + W(y) on K2, where W is a standard Brownian motion. 
We use the notation of section 12.21 Let To be the germ of triangulation of K2 
consisting of the natural simplex To = T(0,e, 1) and Tq its symmetric with 
respect to (1/2, 1/2). In particular, % is a triangulation of K2 itself. Relation 
(fT4|) is satisfied with the network I?. According to proposition 11.11 and theorem 

-^£Xi r ^ -> GL 2 . 

Now we will see that with another germ of triangulation that is not a triangula- 
tion of K2, which in particular does not admit Z 2 as a network, the asymptotic 
convex rearrangement is different. 

Call p the clockwise rotation of M 2 with angle ^ , and consider the new germ 
of triangulation T — p(To) defined by 

' U= (75(^1 + e 2),^(-e!+ e2 )), 
< T = T((0,0),u,l), 
T'=T((0,V2),-u,l), 

r = {T,r}. 

The set of triangles T fulfils condition (|T4|) . and admits Tj- = u\l + M2Z as 
network (as well as \f2e\L + u{L 1 and others). 

Theorem 3.3. Put b 2 n = y/2n. Then, 

pl % =>■ p almost surely, 
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where fi is the centred Gaussian distribution with covariance matrix 
in basis u. So we have, according to theorem ] 1 . 6 A 



2 1 
1 2 



1 

— — 

O n b n 



a.s, —WlVX n -> M,,, — €X n ^C„, 



on int(Kd). 



Proof. Let us take n £ N and z = (x, y) <E Let us write the coordinates of 

(T) 

VXA in the basis u: 

= W + ^) + W (» + ^) - W(x) - W(y), 
X n T i(z) = w(x- -=r-) + W (y + - W(x) - W(y). 



For T': 



= W (x - -L) + W (y- -Lj - - W(y). 



-(vxTw)-^(; J), *.-v>£ 



Hence we set 

62 = V2n, 



A(z) = ^ ^ 



Since cov(y n (z)) converges to A(z) almost everywhere on X^, the candidate /i 
for the limit distribution is the centred Gaussian distribution with covariance 
matrix A(0) in basis u. We now have to verify hypotheses of corollary 12.11 to 
have the convergence /i„ =>■ fx. 

The sum J2 n W = 5 S n ^ i s finite. The covariance function of X is 

t/), (x' ,y')) = xAx'+xAy' + yAx' + yAy'. 
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Using notations of corollary 12.11 we have 

= {(x, y), (x', y')\x = x' or x = y' or x' = y or x = y'}, 

|eni/,f| = o(n 3 ). 

Concerning condition [201 we have to majorisc ||cov (v X„ S \x, y)j | for S G 

{T, T"} when x and y are close. Instead of studying all the cases, we make 
the following remark. This is the covariance matrix of the variables <5i + 62 and 
63 + 64, where the 5{ are Gaussian variables with variance Hence, whatever 
is the dcpcndancc between <$j's, in all cases, each component of the covariance 
matrix will be bounded by — which is the uniform bound we were looking 
for. So, using corollary 12. 11 

a.s, => ,u, 

and the result is proved. □ 

This time, the triangulation matters in the limit distribution. Seeing the 
field under another angle, by rotating the simplexes that generate the trian- 
gulation, induce a radical change in the gradient distribution of the polygonal 
regularisation. More precisely, if p is the rotation applied to triangulation 7o to 
get triangulation T, it does not suffice anymore to rotate p, under p to have the 
new limit measure. This phenomenon illustrates the complex nature of rear- 
rangement for multi-dimensional fields, in contrast with one-dimensional fields. 



3.3 Chentsov field 

This section is devoted to the study of Chentsov field. As the Levy field, it is 
an irregular centred Gaussian field. Like in the previous example, the nature 
of the asymptotic convex rearrangement strongly depends on the choice of the 
triangulation. For z and £ two elements of M. d , denote by z A ( the vector 
which coordinates are the point-wise minimum coordinates of z and Ci an d z 
is the product of coordinates of z. The Chentsov field is defined on as 
the Gaussian field with covariance function F(z, Q) = z A £. Here, we use the 
notation of section 12.21 where X is a Gaussian field. 



Theorem 3.4. Let u = (u l )f =1 be an orthonormal basis and T = T(0,u, 1) a 
regular simplex of (R+) d . We define 

u id = u l A u j - u' A - u 3 A £ R d , i, j £ {1, d} 7 

l(z) = (z 2 • • • Zd, ziz 3 ...z d ,...,Zi... Zd-i), z g K d , (24) 

A^\z) i , j = (l(z),u i , j ), i,je{l,...,d}. 

We call /i( u ) the mixture of the Gaussian probability measures with covariances 
A( u )(z), z g K d , and ip( n > its characteristic function. We have 
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Proof. We will use proposition ^. 21 to compute the only possible limit and corol- 
lary [2J] to show the almost sure convergence. 



Let z £ H^ T \ and l<i,j<d. Then 



E(VX„. 4 (z)VX„.,(z)) = n 2 E((X(z + -u l ) - X(z))(X(z + -xf) - X{z)) 

n n 

= n 2 ( (z + -u*) A (z + -u j ) -{z+ -u l ) A z - z A (z + -u j ) + z 
\ n n n n 

= n 2 \(z+ -u l A u j ) - (z + -u l A 0) - (z + -u j A 0) + z 

\ n n n 



Consider now the function II on M. d defined by II(z) = z. It admits, for all 
z,h <E M. d , the development 

n(z + h) = Oi + hi)...(z d + h d ) = IL(z) + (l(z),h) + q(z, h), 

where 

l(z) = (z 2 —z d ,z 1 z 3 ...z d ,...z 1 ...z d _ 1 ), 

q(z,h) < l(z)\\h\\ 2 

for some positive continuous function 7. It is hence possible to majorise uni- 
formly 7 by a constant C. Hence, 

E(V1„ i(z)VX n Az)) = n 2 (l(z), -u ; A u j - -u ?; A — -u J AO)) + q (z, 

\ n n n / \ n z 

So, we take 

K = \/n, 
Y n = - — vX n , 



b n 

u l A u 3 -u'AO- u J A 0, 



so that we have 

\\A n (z)-A^(z)\\<^ 
Conditions (fT6|) and ([20| arc satisfied. 

Then, according to proposition 12. 2\ the limit distribution ^ T \, if it exists, has 
characteristic function 

<pW{h) = j e -*< fc ' AW « A >d*. 

JK d 
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Now, we have O = {{z, Q £ Kd \ Zi ^ d,l < i < d}, so that V is of class C 1 on 
(A d ) 2 \6. Since (H n T) Y HO 



= 0(n 1 ), using corollary 12. II we have 



□ 

Let T be a germ of triangulation. Let fjT = (^Ter ■ Then we have, 
with probability 1, 

Finding the expression of C^r is not an easy task, and in general we were not 
able to derive explicit formulas. We present here a tractable expression for the 2- 
dimensional Chentsov field with the germ of triangulation To = {T(0, e, 1), T(0, — e, 1)} 
introduced in section [ 



With the notation of theorem 13.41 we have 





ei,i 


= ei, 




ei,2 


= e 2 ,i = 0, 




e2,2 


= e 2 , 


(- 


e)i,i 


= ei, 


(- 


e)i,2 


= (-e) 2 ,i=0, 


(- 


e)2,2 


= e 2 . 



Notice that in this case, we might as well have taken e as basis of T(0, e, 1) or 
e as basis of T(0, — e, 1), but it does not change the result. We have the relief 
to realise that the result does not depend on the choice of the simplexes bases. 
Hence we set 

l(x,y) = (y,x), 
A«W) - (» I). 

Let C a _b = [— oo, a[x [— oo, b[ be a rectangle of S 2 ,a, b 6 M. We are looking- 
for the expression of the asymptotic convex rearrangement C„r , which is the 
unique convex function null in that satisfies 

I- f e -\(h\ly+h%lx) 

X d ((VC ll T )- 1 (C a . b )) = dh — dxdy, (25) 

Jc a , b Jk 2 Vxy 

= G(a)G(b), 
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where 



G(a) = dh —^dx, a S M. 




It is an increasing bijection from K to [0, 1]. In consequence, we define C^t,, by 



Ci(a;) = / G _1 (t) dt 



o 



C fl r (x,y) = (C 1 (x),C 1 (y)), 



Since Ci is convex, so is C„r . We have 



fJ-VC r (C a ,b) = / 'i-c' 1 {x)<a'i-c' 1 (y)<bdxdy = / l G -i (;r .)< a l G -i( y )< b da;(iy 




= G(a)G(b). 



C^tq indeed satisfies (|25|) . 

This function is hence the asymptotic convex rearrangement of X. As in 
example 13. 21 . it can be represented as the sum of two versions of a 1-dimensional 
function 



Nevertheless X cannot be written as the sum of two processes X = Xi X 2 , 
each having asymptotic convex rearrangement C l5 on the contrary of the ex- 
ample at section 13.21 Also, this (J)-factorisation is the exclusivity of the 2- 
dimensional case, and is due to the fact that l(z), defined in (|2~i|) . does not 
involve products of coordinates ZiZj if and only if d < 2. 

Also we show that the asymptotic convex rearrangement of X is in some sense 
the (J) —sum of the average convex rearrangement of the integrated Brownian 
motion. Indeed, let X\ be an integrated Brownian motion. It does not have a 
deterministic convex rearrangement since it is a smooth random function, but 
its gradient distribution satisfies 



E(i v(Xl(BXl) (B) = Mv(Cl0Cl) (B) = ^vc M r (C 0l6 )(5), BeB d . (27) 



Proof. With the previous notations, for a cylinder B = [— oo, a[x [— oo, 6[, since 
X[ is a standard Brownian motion, 



C M r = Ci © Ci. 



(26) 




= G(a)G(b), 



which concludes the proof. 



□ 



23 



The one-dimensional vertical and horizontal restrictions of 2-dimensional 
Chentsov field X are re-scaled Brownian motions, but due to the statistical de- 
pcndancc of distant points, its asymptotic convex rearrangement looks more like 
that of a sum of one-dimensional integrated Brownian motions, in the sense of 
P7|) . Hence the long-range dependancy effects can also have a strong influence 
on the asymptotic convex rearrangement. 



4 Discussion 

In this article we developed tools for computing the asymptotic convex rear- 
rangements of some random fields. We observed that there was a strong de- 
pendancy on the choice of the triangulation used for approximating the field. 
In [5], it becomes apparent that for some 1-dimensional Gaussian processes, 
the Lorenz curve seems a "universal" asymptotic convex rearrangement, in the 
sense that it is the same for polygonal and convoluted approximations . In the 
multivariate case, one would like to do the same, associate to each field a unique 
convex object. For the Levy field and the Gaussian field studied in section [3T21 
(a 0-sum of Gaussian processes) , the range of possible limit distributions seem 
to be a wide class of nondegenrate Gaussian distributions, and it would be in- 
teresting to find a structure in it. This is an ambitious program, and the next 
step is to investigate asymptotic convex rearrangement with other methods of 
approximation, such as convolution, and other ways of representing the asymp- 
totic convex object, like the associated convex body, or the blaschke sum of 
small convex sets involved in the graph of the field approximation. 



5 Proofs 

5.1 Proof of theorem 11.21 

Without loss of generality, we suppose /„ and / convex. It allows us to omit 9Jt 
and £ in the writing. 
((6]) clearly implies ([7]). 

The implication © © follows from the following lemma. 

Lemma 5.1. Let K be a compact convex set, and (/„) a sequence of convex 
functions that converge point-wise to a convex continuous function f on K. 
Then V/ n converges to V/ for the L 1 norm on each convex compact subset of 
int{K). 

Proof of lemma \K7\ We will prove the lemma in three steps. 

Equilipschitzian convex functions on [0, 1]: For k > 0, we put C K the set 

of ft-lipschitzian convex functions on [0, 1]. We are going to show the result in 
the case where / and the (/„) are in C K . In this case, we pick a dense countable 
subset S = {xk 7 k e N} in [0,1]. Since the f' n are bounded (by «), by the 
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diagonal sub-sequence method, we can find a sub-sequence f' a f^ such that, for 
all k, f'a^ixk) converges to some value g(xk), where g is increasing on S. Call 
also g its unique right-continuous increasing continuation on [0, 1], we will show 
that f' n converges to g on [0, 1] for the C 1 norm. Let x be a continuity point of g 
and e > 0. Then, let y and z be in S such that \g(y) — g{z)\ < e and y < x < z. 
For n large enough, since the f' n are increasing, 

-2e < g(x) - g(z) + g(z) - f n {z) < g(x) - fjx) 
<g{x)-g{y)+g{y)-f n (y)<2e. 

Hence f' n converges to g in each of its continuity points, i.e almost everywhere 
according to Riesz-Nagy theorem. Since g is bounded (by re), f' n converges to g 
for the C 1 norm (Lcbcsgue theorem). By integration, g equals /' a.e and so we 
have the result. 



Convex functions on [0,1]: We drop here the assumption that the /„ are 
equilipschitzian. Let / = [a, b] be a compact subinterval of ]0, 1[. Then, for each 
/„, for any x in /, we have, by convexity, 

/n(o)-/n(0) . ^ fn(l)-fn(b) 

1 ^ /«(*) ^ — YZi— 

Since the left and right hand terms converge to finite values when n goes to oo, 
the /„ are equilipschitzian on /, and using the previous result, f n converges to 
/' for the C 1 norm on /. 

Convex functions on K: Let Ij,l < j < d be compact intervals of R 
such that C = 1\ x ... x 1^ is a compact rectangle contained in \x\t(K). Let 
us put, for z in I\ x ... x Id-i, I z — {z} x 1^. For such z, 1 < i < d, define 
G n ,z,i(x) = V/ n ,i(z,x) - Vfi(z,x) on I d , and C n ^(z) = HGn^^H^ . Now we 
have, with Fubini's theorem, 

l|V/„-V/||£i =J2 / C nti (z)dz. 

i=l J hx---xld-i 

Since /„ uniformly converges to / on K, it also does on a segment J z which 
interior contains I z . Each integrand C n ^z) converges point-wise to due to 
the previous result. To dominate it, we write I z =: [a z ,b z ], and call c z a point 
in I z where the monotone function V/ n| ,(z, ■) reaches 0, or c z — a z (arbitrarily) 
if is not reached. Then, using the monotonicity of V/ ni i(z, •) , we have 



C nti (z) < \\fn(a z )\\ + \\f n (b z )\\ + 2\\f n (c z )\\ < 4\\f n \\£ < 4||/||£ + o(l). 

The last upper bound is due to the fact that the point-wise convergence of f n 
to / on the convex C yields uniform convergence. So, Lebesgue's theorem gives 
us the conclusion: 
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||V/ n - V/H&-X). 

Now, each convex compact subset C of int(i^) is contained in a finite union of 
such rectangles, and we have the conclusion. □ 

proof of © ©. 

Let (K e , e > 0) be an increasing family of compact convex sets which satisfy 

|J K e = mt(K). 

e>0 

Let A be a continuity set for /j, and let a > 0. Take e > such that fi{Kf) < a, 
and n after which \fj,^ft(A) — /i^/fe(A)\ < a (A is also a continuity set for /iy/< 
because fiyf^(dA) < fiyf(dA) = 0.). Then, 

|MV/„(A) - /i V /(A)| < A d (A' e c ) + |/iv^(A) - Mv/.(^)| < 2a, 
and it shows that fi n =>■ /i. 

Let us now show ([5]) =4> ((6|). 
This result comes from the structure of convex functions, and of their gradient, 
the "monotone functions" , so we first state a result that helps us apprehend the 
topography of a monotone function. 

Lemma 5.2. There is a class (K e ) e> o of closed subsets of K, satisfying 
(i) e>e' ^K e C K e> , 
(H) U e >o^ = mt{K), 

(Hi) For any convex function f , positive number A and e > , 
P|| V /||([A,ooD < e => Vz g A' e , ||V/(s)|| < 2A 

Hence one can control the locations of points where /'s gradient reaches high 
values. In particular, ||V/|| cannot be "too large" far from the edges of K. 

Proof. Any convex function / on K satisfies 

Vz,CeK, (v/(z)-v/(0,z-0 >o. 

It readily follows from the fact that the restriction of / to [z, (] is convex. Now, 
for z G K, u G K d we introduce the affine cone 

Z(z,u) = {y e K d | (y-z,u}> -\\z - j/||||w||}. 
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Wc have the property that 

y€Z(z,Vf(z))^\\Vf(y)\\>±\\Vf(z)\\. 
Indeed, let y be in Z(z,Vf(z)). 
I|V/(j/)||||j/-z|| > (Vf(y),y-z) > <yf(z),y-z) > ±\\Vf(z)\\\\y- z\\. 

That means that the y in the cone Z{z, V/(z)) cannot have a too small 
gradient, due to the monotonicity property. 

Now we set e(z) = mi u€S d-i Xd(Z(z, u)), which simply plays the role of a lower 
bound for X d (Z(z,V f(z))). We have, for z € K, 

U{y G K | ||V/(y)|| > l -\\Vf{z)\\}) > \ d {Z{z,Vf{z)) > e{z). (28) 
Now we set, for e > 0, K e ={z£ K \ t{z) > e}. For z in K e , 
\\Vf(z)\\>2A^ mwn ([A,oo[)>e(z). 
Hence, given any positive number A, if V/ satisfies 

H\\Vf n \\([A<x>[) < a, 
then, according to (|2^|) . it follows that for z £ K a 

e(z) > a, and so ||V/(z)|| < 2A. 

□ 

We hence have to show that /„ converges to / on int(-ftf). In a first time we 
will use Ascoli-Arzcla theorem to show that the f n uniformly converge on every 
K e , and by consistency they converge point- wise on int(-fT). Then we will show 
that the limit can be nothing but /. 

Since ^v/„ weakly converges to [if, it is a tight family of measures. For all 
e > 0, wc can find A > such that, for all n in N, 

M||V/„||([A<»D < e- 

Hence, according to lemma 

Mn e N,Vz e K e , ||V/„(z)|| < 2A. 

According to Ascoli-Arzela criterion, we know that for all e > 0, {f n (z); z £ K e } 
is a relatively compact family for the uniform convergence. Now, let e be a 
positive number. There exists a convex function f e and a sub-sequence / v s( n ) 
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such that f v t{n) —> fe uniformly on K e . We will show that f e coincides with /, 
which means that / is in fact the limit as only possible limit for a sub-sequence. 
Set efc = j,k > 1 . We build by recurrence <\>\ in the following way. Initiate 
by 4>\ = ¥> c . For k > 1, since (f^ ( n ))n>i is relatively compact, we can find 
<p% +1 such that o V « +1 (n) uniformly converges to a limit fl +l on K ek+1 . The 
/| satisfy the consistency property 

VA; < k' , fl and f%, coincide on K tk . 

Hence we set, for z e Uk>iK ek = mt(K), f e (z) = f%.{z), which does not depend 
on the k such that z is in K ek due to the consistency property. 
Now we set 4> e (k) = <j) e k (k). The sub-sequence (j) e satisfies 

Vzeint(A'), / . (fc) (z)-> f(z). 

Note that we also have the uniform convergence of f^(k) to / e on K e , and so f e 
and / e coincide on K e . By consistency we can build / on int(/\ ) such that /„ 
converges point-wise to / on mt(-fT). According to the result ((6]) =>■ ([5]) proved 
earlier, we know that /iv/„ =^ My/' ana - so ' unicity of the limit, fi v j = /iy/- 
Hence V/ and V/ are two monotone functions on K whose distributions coin- 
cide. The uniqueness in Brenier's theorem ensures us that they are equal a.e. 
We have proved that any cluster point / e of (f n (z),z £ K e ) is equal to / on 
K t . Hence / is the limit of /„ for the uniform convergence on K t . Since for 
convex functions on a convex compact set, uniform convergence and point-wise 
convergence are equivalent, we have 

Ve>0,||/„(z)-/(*)||£^0 

which yields the result. 

5 . 2 Proof of theorem EH 

For the sake of clarity, we omit exponent (T) in this proof. (This concerns 
variables VX„, Y„, fin,ip n )- 
Take h in R d . 

Hereafter, the notation Z refers to a quadruple Z = (z 1 , z 2 , z 3 , z 4 ) <G (Kd) A - We 
set E\ = {1, 2, 3, 4}, V is the class of all subsets of E4, and for < c < 4, V c = 
{P e V I \P\ = c}. We also put £l = e 2 = 1, £3 = e 4 = -1. 



E(\ip n (h) - E(^„(^))| 4 ) = EM <p n (h) - E(tp n (h)j) - E((^„(/i)) 
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E[ [ dz 1 [ dz 2 [ dz 3 [ dz 4 f[ ( e <h,e q Y n (z"))_ e -l(h,cov(Y n ( Z "))h) 
\JK d JK d JK d JK d q=x 

[ E ( y^(_i)l^l e t E«6i»< h ^ y *(**)>-iS s /^(h,oov(y„(««'))h) , \ dZ 

= I V(_i)i p i e ^K' l '( cov (5: 5e p e ^"( z? ))+s: g '?p cov ( y " ( ^' ) ))' l >dz 



'(KaV Pev 



Let us make more precise the term cov (^2 q€P e q Y n (z q ) s j . Let i,j be in 
{l,...,d}. 



E eP ^( zq )) . . = E E t* Y nAz q ) E 

\g£P <?'GP 



gGP q.q'eP 
q^q' 

Thus we define, for P in P, Z in (i^d) 4 , i,j in {1, d}, 

q^q' 

We put temporarily, for Q a finite subset of K d , = -?(K<x*r<y n {z))h) _ 

Notice that for two disjoints subsets Q and Q' of iQ, ipn V'n = V'n U< ^ • Also, 
for Z = {z l5 1 < i < 4} in (Kd) 4 and P C P 4 , we put Z P = {z« | <? 6 P} and 

write for short $f = V^ 4 = ^l* 1 '* 2 '* 3 '* 4 . We have 



E(|^„(/i)-E(93„(/i))| 



E (-l) |P| ^f^ c e^ ^'^-'dZ, 

{K d y Pev 



[ £(-i) |J,| i^(i + E ^xfAzk, + o( x ?(z)r))dz, 

•>{Kd) i PeV l<i,j<d 
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/ ^(E(- 1 ) |P| + EH) |P| E h^wiA dz 

J ( K d) i \PeV PeV l<i,j<d J 



+ 0\ ttxZWdZ • (29) 

\J(K d r J 

Remark that 

E(-d |p| - EE(-D 9 

Per 9=0 Pev q 

= 1-4 + 6-4 + 1 = 0. 

The non-negligible term is a discrete integral, that we can write in the form 
H^J . We put t n = X d (jT) = ^a-Ad(T). For any i, j 

in {1, ...,d}, 



"'(A'd) 4 peV z e (JfT)4 Pe -p 



where \c n \ < Cx/n a.s, according to (fl~5j) . Since Xn(Z)i,j is a sum of terms of 
the form 'E(Y ny i(z)Y n j(C i )) with z,( € the idea is to count the number of 
times each one of these terms appears. 



30 



E ^£(-i) |p| x^k, 

= E(- 1 ) |F| Ew E ^(^(z^^'mi 



PGP q,g ; e-P Ze(if£) 4 



= E(- 1 )' P| E W E E l z ^^(YnAz)YnAO)*Pn 



E(-i) |p| E ^ E v(Y nti (z)Y nJ (o e i,.=,,^=c^f 



PGP ^'eP z,<£H% Ze(HT) 4 

q^q' 

= e E^W^lOEi- 1 ) 1 ' 1 E E ^n' C ' z ' X ' 

q^q' 

I \ 

E ^ } E(Yn,i(z)YnAO) E ^ E (- 1 )' P| 
Vz.CeflT i<,,,'<4 P£P / 

' E ^ z '' c ' } 

Let q q' E E4. Then there are exactly one P of P2, 2 sets P in P3 and 1 
set of P4 that contain q and g'. Hence 

E (-i) |p| = 1-2+1 = 

P<=lT 
P3l,<t' 

and the previous sum is null. The first order term of (|29[) reduces to the deter- 
ministic vanishing non-random term c„. 
Let us estimate the last term of ([29)) . 

i> z nX p n {Z)\\Z 

f E E ^e q e q ,e p e p/ E(Y n 4z^Y n Az q '))E(Y n 4znY n Az p '))dZ, 

= o[ [ 4>n sup|E(r„. 4 (z 1 )y nj (z 2 ))|su P |E(y Il , l (z 3 )y nj (z 4 )|dz ] , 
= o[[/ ^ Z,C} sup |E(y IM (z)y nj (c)|dzdC 

' (Kd) 2 hi J 

And hence formula (fT7| is proved. 
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